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Abstract. We present an extraction of the Sivers distribution functions from the most recent exper- 
imental data of the HERMES and COMPASS experiments, assuming a neghgible contribution of 
sea quark Sivers functions. 
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^ ' In Ref. [1], we presented an extraction of the Sivers distribution functions based 

^ ■ on a fit of SIDIS experimental data from the HERMES [2] and COMR\SS [3, 4] 

^ ■ collaborations. Data from HERMES [2] presented an unexpectedly large A^Jf ' 

— j. . asymmetry for K^ production, about twice as much as the analogous asymmetry for 

1^ I n^. Such a large asymmetry suggested an important role of the sea Sivers functions. Our 

O ' analysis confirmed this expectation finding a large contribution of the 5-Sivers function. 

Since then, new experimental results have become available: the COMR^SS collab- 
oration has released new SIDIS data off a proton target, showing a clear Sivers asym- 
metry [5]; a new HERMES data analysis, based on a much larger statistics, while con- 
^ ■ firming the previous pion data, softens the enhanced peak in the K^ Sivers azimuthal 

moment [6]. These new data prompted us to perform a new analysis. Here we present a 
preliminary phenomenological fit of COMR^SS [3, 4] and HERMES [6] data including 
only u and d valence quarks in order to see if the new data can be described consistently 
without any contribution of the sea Sivers functions. 

The SIDIS transverse single spin asymmetry (SSA) Ayj measured by HER- 

MES and COMR\SS, in the f - p cm. frame and at order k^/Q, is given by [7, 8, 9]: 

sm{^,-^s) _ ^qI^^sd(l)hd^k^A^f^/pt{x,k^) sin((p - ^s)^^^D\{z,Pa.) sm{(j)h - (^5) 

(1) 
where ^s ^nd 0/, are the azimuthal angles identifying the directions of the proton spin S 
and of the momentum of the outgoing hadron h respectively w.r.t. the lepton plane, while 
cp defines the direction of the incoming (and outgoing) quark transverse momentum. 



k^ = k^{cos(p,sm^,0); '^^ ^ is the unpolarized cross section for the elementary 

scattering iq — > iq; Dq(z,p±) is the fragmentation function describing the hadronization 
of the final quark q into the detected hadron h with momentum F/,; h carries, with respect 
to the fragmenting quark, a light-cone momentum fraction z and a transverse momentum 



p^. Finally A^f , i^{xjk±) is the Sivers distribution function [10]: 



2k \ 

A%/p,{x,k^) = -^/iV(^,fc^) . (2) 

nip 

The Sivers function is parameterized in terms of the unpolarized distribution function, 
as in Ref. [7], in the following factorized form: 

A%/pt{x,k^) = 2^q{x)h{k^)f^ip{x,k^) , (3) 

with 



:NqX^^{\-xf^^-^^^-^ , and h{k^) = Vie 



^a{x)=Nax'^H\ -xf'i ' H ' r-HJ ^^^ hikA = V2e^e 



-k\/M- 



<'i8,^^ ' ^1 



(4) 

where A^^, Uq, j8^ and M\ (GeV/c) are free parameters to be determined by fitting 
the experimental data. We adopt a Gaussian factorization for the unpolarized distri- 
bution and fragmentation functions with the Gaussian widths (fc^) and (p^) fixed 
to the values found in Ref. [8] by analysing the Cahn effect in unpolarized SIDIS: 
{k\) = 0.25 (GeV/c)2 and {p\) = 0.20 (GeV/c)^. For the unpolarized, A:j^ -integrated 
distribution and fragmentation functions we use the GRV98 [11] and DSS [12] sets. 
We best fit the HERMES proton and COMPASS deuteron data from Refs. [3, 6] (209 
points) including only Sivers functions for u and d quarks, corresponding to seven free 
parameters, shown in table 1. The results we obtain are rather satisfactory, with a jj^^ 
of about 1.06. They are shown in Figs. 1, 2. The corresponding Sivers functions are 
plotted in the right panel of Fig. 3. Such results are similar to those obtained in Ref. [7] 
The gray band in Figs. 1, 2, 3 represents the statistical error of the fitting procedure, 
calculated as in Ref. [1]. As shown in the left panels of Figs. 1 and 2, pions data are 
well described by our fit. The new HERMES data on kaon production can be described 
reasonably well without any sea Sivers functions contribution, although their inclusion 
in the fit procedure can slightly improve the Xdof [^^]- ^'^ ^^^ ^^^^ panel of Fig. 3 we show 
our predictions for the Sivers asymmetry at COMR^SS kinematics on a proton target. 
Although apparently we overestimate the positive charged hadron asymmetry, COM- 
R^.SS proton data are affected by a scale error of ±0.01, not shown in the figure, so that 
our predictions are still compatible with them. Figs. 4 and 5 show the uncertainties of 
the fit parameters. They are obtained following the procedure of Ref. [1]. We generate 
randomly new sets of parameters, then we collect the sets giving a variation of x^ with 
respect to the minimum corresponding to a 95.45% confidence level. The figures show 
the Ax^ as a function of each parameter. Some of the parameters are correlated. The 
strongest correlation is between Nu and M^, as it is shown in the most right panel of 
Fig. 5. It is interesting to notice that the Mi parameter, which fixes the Gaussian width 
of the Sivers function (i.e. its distribution in k^ is rather well constrained and turns out 
to be between one half and two thirds of the unpolarized distribution function width. 
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FIGURE 1. Fit of HERMES data [6] for pion (left panel) and kaon production (right panel). 
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FIGURE 2. Fit of COMPASS deuteron data [3] for pion (left panel) and kaon production (right panel). 
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FIGURE 3. In the left panel we compare our predictions on proton target for charged hadrons with the 
data released by the COMPASS collaboration [5]. The errors on the data are the statical and systematic 
errors added in quadrature. The right panel shows the first moment of the Sivers functions extracted from 
the fitting procedure. 



TABLE 1. X and best values of the parame- 
ters. 

x^ldof^im 

Nu = 0.40 a„ = 0.35 A- = 2.6 

Nd = -0.97 ad = 0.44 jS,/ == 0.90 

M2=0.19GeV2 



68% C.I. 
95.45% C.I. 
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FIGURE 4. From left to right: A^^ as a function of A^^,, Nd, cCu and a^/. 
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FIGURE 5. From left to right: A;^^ as a function of j3„, j3^/, M\ and correlation between M\ and A^^,. 



